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The magnetic and thermodynamic properties of a new ferrimagnetic decorated spin-(^, 1) Heisen- 
berg chain with spin-1 pendant spins are investigated for three cases: (A) J\, J2>0; (B) Ji>0, J2<0; 
and (C) Ji<0, J2>0, where Ji and J2 are the exchange couplings between spins in the chain and 
along the rung, respectively. The low-lying and magnetic properties are explored jointly by the 
real-space renormalization group, spin wave, and density-matrix renormalization group methods, 
while the transfer-matrix renormalization group method is invoked to study the thermodynamics. 
It is found that the magnon spectra consist of a gapless and two gapped branches. Two branches 
in case (C) have intersections. The coupling dependence of low-energy gaps are analyzed. In a 
magnetic field, a m=| (m is the magnetization per unit cell) plateau is observed for case (A), while 
two plateaux at m=i and § are observed for cases (B) and (C). Between the two plateaux in cases 
(B) and (C), the sublattice magnetizations for the spins coupled by ferromagnetic interactions have 
novel decreasing regions with increasing the magnetic field. At finite temperature, the zero-field 
susceptibility temperature product yT and specific heat exhibit distinct exotic features with vary- 
ing the couplings and temperature for different cases. yT is found to converge as which is 
different from the divergent behavior in the spin-(i,l) mixed-spin chain without pendants. The 
observed thermodynamic behaviors are also discussed with the help of their low-lying excitations. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Mg, 75.50.Gg 



I. INTRODUCTION 

In recent years, one-dimensional (ID) quantum fer- 
rimagnets with two kinds of antiferromagnetically 
exchange-coupled centers have attracted much atten- 
tion due to their exotic properties. The large 
families of compounds ACu(pba)(H 2 0)3-nII 2 and 
ACu(pbaOH)(H 2 0) 3 -nH 2 0, where A=Mn, Fe, Co, 
Ni, Zn, pba=l,3-propylenebis, and pbaOH=2-hydroxy- 
1,3-propylenebis, have been extensively explored in 
chemistry,— which are good realizations of the mixed-spin 
chains. These compounds exhibit typically the ID ferri- 
magnetic behavior of yT (x is the magnetic susceptibility 
and T is the temperature) that shows a rounded mini- 
mum with temperature^ 

As a simple model to describe the mixed-spin chains, 
the antiferromagnetically coupled spin-(i,l) Heisenberg 
chain have also been extensively studied by various meth- 
ods, such as the spin wave theory,—^ Schwinger bo- 
son mean field^ density-matrix renormalization group 
(DMRG), 3 quantum Monte Carlo^ and so on^ It has 
been found that its ground state has a spontaneous mag- 
netization at m=^ (m is the magnetization per unit cell) 
that is consistent with the Lieb-Mattis theorem^ and 
the system has a ferrimagnetic (FI) long-range order. 
The one-magnon excitation spectra consist of a gapless 
ferromagnetic (FM) branch from Sg to Sg — 1 (Sg is the 
good quantum number of total spin in z component in 
the ground state) and a gapped antiferromagnetic (AFM) 
branch from Sg to Sg + 1. 11 This magnon gap was numer- 
ically found to be 1.759J (J is the exchange coupling). 3 
In a magnetic field, the system exhibits a magnetization 
plateau at m—^ with the width of 1.759J, corresponding 
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FIG. 1: (Color online) Sketch of the spm-(r, S) decorated 
Heisenberg chain with spin-er pendant spins. 



to the gap of the AFM magnon branch. 12 Different from 
the 5*=! Haldane chain, in this mixed-spin chain the spin 
gap (1.2795J) from the ground state to the lowest state 
in the subspace of Sg + 1 is less than the magnon gap 
(1.759 J) and thus is not a magnon-like excitation. 3 The 
thermodynamic properties in the coexistence of the AFM 
and FM excitations^^ and in the critical phase under a 
magnetic field&i^ have also been investigated. 

Recently, another interesting family of cyanide-bridged 
coordination compounds with pendant magnetic ions 
are synthesized in experiment j 15 ' 16 One of them is 
the cyanide-bridged Ni(II)-Fe(III) complex with a novel 
building block [Fe(l-CH 3 im)(CN) 5 ] 2 -, 17 which can be 
treated as the ID structure as shown schematically in 
Fig. [T] owing to the weak interchain interactions, where 
the Ni(II) (Si and Oj) and Fe(III) (r«) ions have spin 1 
and 4 , respectively. This compound realizes a decorated 
spin-(i,l) mixed-spin chain with spin-1 pendant spins. 
Although the intrachain couplings Ji<0 and J 2 <0 are 
both FM interactions in the present compound, it is no- 
ticed that any other couplings (i.e., Ji, J 2 >0, Ji>0 and 
J 2 <0, Ji<0 and J 2 >0) would give rise to ferrimagnets, 
making the realization of such a FI structure more ac- 
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cessible to the experiment. This family of mixed-spin 
chains with pendant spins provides a new scheme to 
study the ID quantum ferrimagnetism, which may have 
new exotic properties. Although the influences of pen- 
dant spins on magnetism have been discussed in some 
antiferromagnets,— the studies on such ferrimagnets are 
still rare till now. 

In this paper, we shall explore the physical proper- 
ties of this new ferrimagnetic structure, and compare 
with the spin-(i, 1) mixed-spin chain without pendants. 
The low-lying, magnetic and thermodynamic properties 
of the spin-(i, 1) decorated Heisenberg chain with spin-1 
pendant spins for three cases: (A) Ji, J 2 >0; (B) Ji>0, 
J2<0; and (C) Ji<0, J 2 >0 will be studied using vari- 
ous techniques. It is unveiled that due to the pendant 
spins, the new ferrimagnets exhibit rather distinct mag- 
netic and thermodynamic behaviors from those of the 
spin-(i, 1) mixed-spin chain. The three cases with differ- 
ent couplings are uncovered to have their own features, 
which are expected to observe in experiments. The exotic 
properties of this system will also shed light on further 
understandings of quantum ferrimagnetism. 

This paper is organized as follows. In Sec. II, the 
model Hamiltonian is introduced. In Sec. Ill, the low- 
energy effective Hamiltonians in both strong and weak 
couplings are analyzed utilizing the real-space renormal- 
ization group (RSRG) method. The low-lying excita- 
tions and magnetic properties are investigated by the 
linear spin wave (LSW) theory and DMRG in Sec. IV. 
In Sec. V, we shall study the zero-field thermodynamics 
by means of the transfer-matrix renormalization group 
(TMRG). Finally, a summary and discussion will be 
given in Sec. VI. 

II. MODEL HAMILTONIAN 

The Hamiltonian of the spin-(i,l) decorated Heisen- 
berg chain with spin-1 pendant spins in a magnetic field 
can be written as 

N 

H = Jin ■ Si + J\S l ■ fj+i + J 2 n ■ di) 

i=l 
N 

- hJ2(T*+S* + a*), (1) 
t=i 

where Si is the a=l pendant spin, ?i and o, are the spins 
in the chain with r=i and 5=1, respectively, Ji,2>0 
(<0) denote the AFM (FM) couplings, and h is the mag- 
netic field. Throughout the context, we take J\ as an 
energy scale and gfiB=^- The schematic representation 
of the model is shown in Fig. Q] 

Analogous to the spin-(i,l) mixed-spin chain with- 
out pendants, the system with Hamiltonian ([1} has a 
spontaneous magnetization in the absence of magnetic 
field according to the Lieb-Mattis theorem^ In case (A) 
(Ji.2>0) the spontaneous magnetization per unit cell is 



m=|, while in both cases (B) (Ji>0, J2<0) and (C) 
(Ji<0, </2>0) it is spontaneously magnetized at m=\. 
The Goldstone theorem—^ allows gapless excitations in 
these cases owing to the spontaneous breaking of the 
SU(2) symmetry. 

III. REAL SPACE RENORMALIZATION 
GROUP ANALYSIS 

In this section, the low-energy effective Hamiltonians 
of the three cases in both strong- and weak-coupling lim- 
its are derived utilizing the RSRG. 20 In the RSRG pro- 
cedure, the Hamiltonian is divided into intrablock (H B ) 
and interblock (H BB ) parts. By diagonalizing H B , a 
number of low-energy states are kept to project the full 
Hamiltonian into the renormalized Hilbert space. Al- 
though RSRG cannot give the results as accurate as the 
numerical approaches, it can give a good qualitative de- 
scription for low-energy properties. 

A. Ji>0, J 2 >0 

Let us first consider the strong-coupling limit ( J2^> Ji)- 
Since the interaction between n and cr.; is strong, each 
rung can be considered as the isolated block in the first 
step of renormalization group (RG). Each block consists 
of two multiplets whose total spins are 1/2 and 3/2 with 
energies — J 2 and J2/2, respectively. The spin-i doublets 
are kept as the basis to construct the embedding opera- 
tor T to project the full Hamiltonian onto the truncated 
Hilbert space. The effective Hamiltonian can be obtained 
as 

1 N 

K T ° m ' = -NJi ~ ^1 £(3 -Si + Si- S' i+1 ), (2) 

6 i=l 

where 5^=1/2 is the renormalized spin truncated from 
the rung block. The Hamiltonian © describes a spin- 
(|,1) mixed-spin chain with a renormalized FM coupling 
— \J\- In the next step of RG procedure, the Hamilto- 
nian (0) is further projected onto a 5"=3/2 FM Heisen- 
berg chain 

T 9 N 

iJ« = _ ( j 2 + ^ )N ±_ Jx £ s>> ■ S? +1 . (3) 

i—l 

The magnon excitations in this FM Hamiltonian corre- 
spond to the magnons from Sq to Sg — 1 in the original 
system, which are hence expected to be gapless with a 
quadratic dispersion in low energies. The RG can also 
give the sublattice magnetization ms=l, m r =—^, and 
m CT =|. The sum of them gives |, recovering the sponta- 
neous magnetization of the original system. 

In the weak-coupling limit ( J 2 <Ct/i), the spins and Si 
are taken as a block, and the doublets with spin 1 /2 are 
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kept to truncate the block Hilbert space. The effective 
Hamiltonian is 

N 

£wca k ,l = _ NJi _ ^(-Jj' ■ S' i+1 + -J 2 S' t ■ Si), (4) 
i — 1 

where S[= 1/2 is the renormalized block spin. The Hamil- 
tonian is mapped onto a spin-1/2 FM Heisenberg chain 
with a—1 pendants coupled by the renormalized FM in- 
teraction — |Ja- The spin wave analysis unveils that it 
has a gapless FM excitation with the dispersion 

^fe ~ ^Jik 2 (5) 

for fc— >0, which corresponds to the magnon excitation 
from Sq to Sq — 1 of the original system. The sublat- 
tice magnetization is obtained as ms=§, m T =—^, and 
m CT =l, and the sum of them is also §. 

From Eqs. © and ([5]) it can be seen that in the two 
coupling limits, the low-energy behaviors of the gapless 
branch are both dominated by J\. Besides, ras and mo- 
exchange their values in the two limits, and m T is un- 
changed, implying a possible crossing of m T in the inter- 
mediate region of J2/J1, which would be confirmed by 
the DMRG results in the next section. 



B. Ji>0, J 2 <0 

In the strong-coupling limit (IJ2 I^^Oj owing to the 
strong FM J 2 , the low-energy multiplets with total spin 
3/2 are kept to project the Hamiltonian in the first step 
of RG. The effective Hamiltonian is obtained as 

1 1 N 

^strong,! = 1 + ^ ^ . g , + g , . ^ (g) 

1—1 

with the renormalized spin of rung block ^'=3/2, which 
depicts a spin-(|, 1) Heisenberg chain with a renormal- 
ized AFM coupling \J\. In the next step of RG proce- 
dure the Hamiltonian ((6]) is projected to a S"—l/2 FM 
Heisenberg chain 

1 m N 

Hf^ = i( 3 J 2 - 5J x )iV - g J X £ £/ • S^, (7) 

i— 1 

whose FM excitations imply that the magnon excitations 
of the original system from Sq to Sq — 1 are gapless 
with a quadratic dispersion relation in low energies. The 
sublattice magnetization is obtained as ms=-|, m T = j^, 
and to ct =|, giving the spontaneous magnetization per 
unit cell m=-i. 

In the weak-coupling limit (|</ 2 |<SJi), we perform the 
RG on the block spins Si and Tj with the doublet of total 
spin 1/2. The effective Hamiltonian is 

N 

^wcaM = _ NJi _ YQjJ' ■ S' l+1 + l -J 2 S> ■ Si) (8) 



with the renormalized spin 5—1/2, which describes a 
spin-1/2 FM Heisenberg chain with a— I pendant spins 
coupled by the renormalized AFM interaction —\J%- The 
spin wave results show that the excitations that corre- 
spond to those from Sq to Sg — 1 in the original system 
are also gapless with 

u k ~\hk 2 (9) 

for fc— >0. In the two limits, it is observed from Eqs. 
(|6]) and ((9|) that the low-energy behaviors of the gapless 
excitation are both dominated by J\. 

C. Ji<0, J 2 >0 

For the strong-coupling limit (J 2 :»|Ji|), because of the 
strong AFM J 2 , Si and fi are renormalized by the doublet 
with total spin 1/2. The effective Hamiltonian is given 
by 

1 N 

^strong,! = ~NJ 2 - -J, ■ £ + £ • S[ +1 ) (10) 

i=l 

with S"=l/2. Eq. (fTU| describes a spin-(i,l) Heisen- 
berg chain coupled by the renormalized AFM interaction 
— 1 Ji. In the second step of RG, the Hamiltonian (ITU)) is 
projected to a S"'=l/2 FM Heisenberg chain 

J A N 

tf« = -(J2 Y )N + 27 Jl ^ § ? ■ § ™> (U) 

i— 1 

which unveils the gapless excitations from Sg to Sg — 1 
of the original system. The sublattice magnetization is 
obtained as mg=|, m T =-^, and m a =— |, giving rise to 
the spontaneous magnetization per unit cell m=^. 

In the weak-coupling limit (J2^C| JiQ, the spins Si and 
fj are taken as a block, and the multiplet with spin 3/2 
are kept to truncate the block Hilbert space. The effec- 
tive Hamiltonian is obtained as 

# wea k) l = 1 ^ + ^ J lS > . S> +1 + l -J 2 S> ■ Si), (12) 
i— 1 

which describes a spin-3/2 FM Heisenberg chain with 
antiferromagnetically coupled pendant spins <ii. The spin 
wave analysis indicates that the effective system has a 
FM gapless excitation with 

uj k ~ -Jifc 2 (13) 

for fc— s>0. Analogous to the above cases, the low-energy 
behavior of the gapless branch in the two limits are de- 
termined by Ji, which can be seen from Eqs. (|10p and 

USD. 

Based on the RSRG analyses, one may observe that the 
different cases have distinct low-energy effective Hamil- 
tonians, and the magnon excitations from Sg to Sg — 1 
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are always FM and gapless, being consistent with those of 
the spin-(i, 1) mixed spin chain. The dispersion relations 
near k=0 are found to be dominated by J\ in both strong 
and weak couplings. The low-energy effective Hamilto- 
nians for cases (B) and (C) in the strong-coupling limit 
are analogous except the magnitude of spin, whose ther- 
modynamic properties will be compared in Sec. V. 



IV. LOW-LYING EXCITATIONS AND 
MAGNETIC PROPERTIES 

In this section, the low-lying excitations and magnetic 
properties are explored by means of the LSW— « and 
DMRG^ During the DMRG calculations, the chain 
length is taken as L=300, and the Hilbert space is trun- 
cated to 240 most relevant states. Open boundary con- 
ditions are adopted, and the truncation error is less than 
10 -8 in all calculations. 



A. Ji>0, J 2 >0 

The Holstein-Primakoff (HP) transformations are in- 
troduced as follows: 



of = si - ajai, 

af = y 2si - a\ai<ii, 

a i = a \\j 2s i ~~ a l a i 
for the sublattice of c?j spins with Si=l, and 

r? = -s 2 + b\b u 



(14) 



r+ = b\J2 S2 -b\b 



2s 2 - blbA 



(15) 



for the sublattice of spins fi with S2=\ 1 where the oper- 
ators ai and bi are bosons. The spins Si are transformed 
in the similar way as Eq. (fill with bosonic operators 
Ci and c\. Thus, the magnon spectra can be obtained 
by diagonalizing the Hamiltonian after performing the 
Fourier and Bogoliubov transformations. As shown in 
Fig. [U(a), the magnon spectra consist of a gapless branch 
u!i t k and two gapped ones u>2,k and u>3 t k- In the presence 
of a magnetic field h, both (Ji t k and u)2,k increase, while 
LU3 t k decreases, indicating that coi.k and w 2i fc describe the 
magnons from Sg to Sq — 1 while u>s t k are those from 
Sq to Sq + 1- For J2—O, the spectra are reduced to a 
gapped and a gapless excitations, which agree exactly 
with those of the spin-(|,l) mixed-spin chain. When J 2 
is set in, the gapless branch splits into uii t k and U)2,k- 
With increasing J 2 , both w 2 ,fe and 0)3^ enhance. It is 
found that the low-energy dispersions near k=0 of the 
gapless branch 0J± t k are insensitive to J2 but dominated 
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FIG. 2: (Color online) (a) Magnon excitation dispersion for 
case (A) with Ji=J 2 =l under different magnetic fields, (b) 
Magnetization curves for different J2. The inset shows the 
coupling dependence of the gap Amp, Asw, and As G + i. (c) 
Local magnetization as a function of lattice site for h<h c i . (d) 
Coupling dependence of sublattice magnetization for J 2 /Ji=l 
and h<h c . 



by J\ in a wide range of the coupling ratio, which covers 
the result obtained from the RSRG. 

The magnetic curve m(h) and low-energy gaps are then 
studied by the DMRG method. As shown in Fig. I^b), 
m(h) has a plateau at the spontaneous magnetization 
m=|, whose width Amp increases with increasing J 2 . 
In this figure, h c \ denotes the field where the plateau 
disappears, and h s is the saturation field. For J2=0, 
Amp reduces to 1.759Ji of the spin-(i,l) mixed-spin 
chain, which is exactly the gap of its massive magnon 
branch i 11 ' 12 The coupling dependence of A MP is illus- 
trated in the inset of Fig. Hfb), showing that Amp in- 
creases almost as a linear behavior. The J\ dependence 
of Amp is also studied by taking J 2 as the energy scale, 
which is not presented here. It is found that Amp/J2 
varies rather slowly with Ji, which means that Amp is 
mainly scaled by J 2 in this case. The gap of the massive 
magnon branch £^3^ (Asw) is also shown in the inset of 
Fig. [H^b) in comparison to Amp- The magnon gap ob- 
tained from the LSW appears to be smaller than Amp, 
where the deviation increases for stronger J 2 . It appears 
that the LSW underestimates the magnon gap from Sg 
to Sq + 1. 

We also compute the spin gap As G +i from the ground 
state to the lowest state in the Sg + 1 subspace, as shown 
in the inset of Fig. Oh). Analogous to the spin-(i,l) 
mixed-spin chain, As G +i is smaller than Amp, indicating 
that Ag G+1 is also not a magnon-like excitation. But, it 
has a similar behavior with coupling to Amp, which can 
thus be used to describe the low-energy behaviors. The 
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FIG. 3: (Color online) (a) Magnon excitation dispersion for 
case (B) with Ji = l and J2=— 1. (b) Magnetization curves 
for different J2. The inset shows the coupling dependence of 
the gap Amp, Asw, and As G +i. (c) Coupling dependence 
of sublattice magnetization, (d) Local magnetization as a 
function of lattice site for h C 2<h<h C 3- 



spin gap from the ground state to the lowest state in 
the Sq — 1 subspace is computed, which is found always 
vanishing and is consistent with the gapless branch u)\ 



Next let us discuss the spin-spin correlation function 
and local magnetization in ground states. Figure (He) 
shows the local magnetization as a function of lattice 
site for h<h c \. In the ground state, the spin correlation 
functions along the chain have a long-range order (LRO), 
and the spin fluctuations (Sf S*)-{S* }{Sj) decay rather 
rapidly (not presented here). Hence it is adequate to 
study only the local magnetization in the ground state. 
For J 2 =0, it gives m T =-0.29248 and rn s =0. 79248, ex- 
actly in agreement with the previous result^ After tuning 
on J2, the coupling dependence of sublattice magnetiza- 
tion is shown in Fig. [2jd). The pendant spin magnetiza- 
tion m a is suppressed by the quantum fluctuations after 
tuning J2, while ras increases and approaches saturation 
for large J2, confirming the results of the RSRG. The in- 
teresting phenomenon is the behavior of m T . As shown 
by the arrow in Fig. [2fd), m T decreases for J 2 /Ji<l, 
and turns to increase when Ji j J\ > 1 , which has a turn- 
ing point at J2/J1—I, as suggested by the result of the 
RSRG. Meanwhile, it can be seen that m a and ms inter- 
sect near J2/ ' J\=l. The changes of coupling dependence 
of the magnetic moments near J<ij J\=\ may be owing to 
the competition of the two AFM interactions. 



B. Ji>0, J 2 <0 

To perform the LSW calculation, 5i and T} spins are 
transformed as the form of Eq. (IT4l) by the bosonic op- 
erators (<2i,<4) and (bi,b\) with si=l and i, respectively, 

while Si are transformed as the form of Eq. (fT5j) by (cj ,cj ) 
with S2—1- As shown in Fig. [3ta), the spectra consist of 
a gapless and two gapped branches. In the presence of 
magnetic field, the gapless branch uii t k and the gapped 
branch increase, while the gapped branch u>2.k de- 
creases, indicating that Wi k and ^3,k are the excitations 
from the sector Sq to Sq — 1 while uj2,k describes the 
excitations from Sq to Sg + 1. It can be seen that the 
spectra in this case are quite different from those of case 
(A). The branch uj2,k from Sg to Sg + 1 is close to the 
gapless branch for the present case, and with increas- 
ing |J 2 |, it increases slightly. The resulting distinctions 
in the thermodynamics would be explored in Sec. V. 
With changing the couplings, it is found that the low- 
energy dispersions of the gapless branch wi fe near k—0 
are dominated by J\ , which is consistent with the RSRG 
result. 

The magnetic curve m(h) and low-energy gaps are 
shown in Fig. [3jb). It can be seen that m{h) exhibits 
two plateaux at m=^ and §. We denote the field where 
the fn^\ plateau vanishes as h&, the lower and up- 
per critical fields for the w=| plateau as h c2 and h c ^, 
respectively, and h s as the saturation field. With in- 
creasing IJ2I, the width of the m=\ plateau Amp is 
enlarged slightly, while the rn=| plateau decreases and 
smears when \J 2 \/ J\~>2.0. It should be noted that in 
some quasi-one dimensional polymerized Heisenberg an- 
tifcrromagnets, there might be a transition from the 
plateau state to the non-plateau state that is usually of 
the Kosterlitz-Thouless type^ Such a transition point 
cannot be numerically determined accurately owing to 
the finite-size length of the chain. Therefore, in the 
present case, whether or not a plateau-non-plateau tran- 
sition with couplings at the m=% plateau exists is still 
uncertain from our DMRG numerical results. The cou- 
pling dependence of Amp is illustrated in the inset of 
Fig. [21b), showing that Amp increases with enhancing 
\J2\1 and different from case (A), Amp goes to saturate 
at large IJ2I, which suggests that Amp is mainly scaled 
by Ji for large \J2\I Ji- The coupling dependence of the 
magnon gap k>2,fc=o (Asw) is also shown in the inset of 
Fig. [21b). It can be seen that the spin wave is capable of 
describing the coupling dependence of magnon gap qual- 
itatively, though it underestimates the value like in the 
case (A). The spin gap As G +i from the ground state to 
the lowest state in the Sq + 1 subspace is also computed. 
As shown in the inset of Fig. [3l As G +i is less than Amp, 
indicating that it is also not a magnon-like excitation, 
and its behavior for different coupling ratios is consistent 
with Amp and Asw- 

The coupling dependence of sublattice magnetization 
in the ground states is shown in Fig. Etc). As the quan- 
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FIG. 4: (Color online) Magnetic field dependence of the sub- 
lattice magnetization ms, m T and m a for |J2|/Ji = l and 2.5. 
The arrows indicate the minimum of m T , where the magneti- 
zation per unit cell is m=3/2. 



turn fluctuations become strong after tuning J2, the pen- 
dant spin magnetization m a decreases with increasing 
J2I, which is analogous to case (A). However, the FM 
coupling has different effects on the spins in the chain 
compared with the AFM J2. With increasing J%\, both 
m T and ms increase slightly, and m T does not show an 
extremum like in the case (A). 

In the m=| plateau region (h C 2<h<hc3), the local 
magnetization is shown in Fig. Eld) for \J2\I Ji—1 as an 
example. It may be expected that all the local magnetic 
moments would increase from h c \ to h c i. However, by 
comparing the local magnetization below h c \ [Fig. GHc)] 
and in the plateau [Fig. [3jd)] , it is surprising to notice 
that m T decreases from 0.3156 below h c i to —0.0046 in 
h c i- Therefore, the field dependence of sublattice mag- 
netization is studied, as shown by ms, m T , and m a for 
\J2\/Ji=l and 2.5 in Fig. |4j It can be seen that m T de- 
creases continuously from h c \ to h C 2 , while to ct decreases 
in a short range above h c \. For a comparison, we also 
calculated the sublattice magnetization as a function of 
field for both the spin-(i, 1) mixed-spin chain and the 
case (A). The results shows that the above behavior is 
not seen. Therefore, this novel decreasing behavior may 
be owing to the competition between the FM and AFM 
interactions in a magnetic field. 



C. Ji<0, J 2 >0 

For Ji<0 and J2>0, the HP transformations are ap- 
plied on the spins Si and fi with the form of Eq. (|T4|) 
by (bi,b\) and (ci,c\) for si—1 and i, respectively, and 
that with the form of Eq. (|15j) is applied on Si for ,S2=1 
with (ai,a\). As shown in Fig. EJa), the spectra consist 
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FIG. 5: (Color online) (a) Magnon excitation dispersion for 
case (C) with Ji=— 1 and Ja=l. (b) Magnetization curves 
for different J2. The inset shows the coupling dependence of 
the gap Amp, Asw, and As G +i. (c) Coupling dependence 
of sublattice magnetization, (d) Local magnetization as a 
function of lattice site for h C 2<h<h c ^. 



of a gapless (wi,fc) and a gapped (u^fe) magnon branches 
from the sector Sg to Sq — 1, as well as a gapped one 
(w2,fc) from 5g to + 1, which can be identified by the 
shifts of the branches with the magnetic field. It is no- 
ticed that the gapped branch uj2,k from Sg to Sq + 1 is 
close to the gapless branch wi^, which is similar to the 
case (B), but W2,k has lower energies than wi^ for large 
wave momenta k in the present case. With increasing 
J2/IJ1I, L02,k enhances and the intersected momenta of 
the two branches shift to higher values. A similar inter- 
section of magnon branches has also been observed in the 
spin-(i,l) mixed-spin chain with AFM nearest- neighbor 
and FM next-nearest-neighbor interactions^ The influ- 
ences of this intersection on the thermodynamics would 
be discussed in the next section. For u>i k, it is found that 
the low-energy dispersions near fc=0 are also dominated 
by Ji, which agrees with the RSRG analysis. 

In Fig. EJb), the magnetic curves m(h) for different 
couplings are shown. Similar to case (B), m(h) has two 
plateaux at m=i and |, whose critical fields are denoted 
by the same symbols as the case (B). With increasing 
^2/|^i|, the width of the m=^ plateau (Amp) extends, 
while that of the m=f plateau is enlarged, which dif- 
fers from the case (B) where the m=| plateau decreases 
with increasing the coupling ratio. The inset of Fig. [5jb) 
shows the coupling dependence of the low-energy gaps 
Amp, Asw, and As G +i- It can be seen that the gaps 
behave similarly to those in the case (B). The gaps ap- 
proach to the saturation for large Jij\ indicating that 
they are mainly scaled by J\ in the large J2 limit. The 



7 



1.0 



0.5 



o 



03 
N 



CD 

c 

« 0.0 

<D 
O 

§ -0.5 

Z3 
CO 



-n-n-H __ D -O^^D-D-D-D-D- 

m=3/2 A 



-1.0 



o-o-o-o-o 
^m=3/2 



O' 



J /L -°-J 2 /|J,|=0.5, m s -n-J 2 /|J l |=1, m s 
-O-J^IJ^O.5, m -o- J 2 /|J,|=1 , m t 



-^■-J/|J,|=0.5, m o -A- ^,1=1, m o 



0.0 



0.2 



0.4 0.6 

h/u I 



0.8 



1.0 



FIG. 6: (Color online) Magnetic field dependence of the sub- 
lattice magnetization ms, m T and m a for J2/|Ji|=0.5 and 1. 
The arrows indicate the magnetic field h C 2 where the magne- 
tization per unit m=3/2. 



LSW also underestimates the magnon gap of ui2,k as Asw 
is smaller than Amp- Ag G +i appears to be smaller than 
Amp, which means that the spin gap from the ground 
state to the lowest state in the subspace with Sq + 1 is 
also not a magnon-like excitation. In the ground states, 
the coupling dependence of the sublattice magnetization 
is displayed in Fig. (5Jc) . It can be seen that as the quan- 
tum fluctuations are induced by J2 , ms and m T decrease 
with increasing J2, while m a increases. In this case, the 
sublattice magnetic moments have more prominent vari- 
ations with the change of the couplings than the previous 
cases. 

In the ?n=! plateau {h C 2<h<h c ^} , the local magnetic 
moments for J2/I Ji|=l are shown in Fig. HJJd). By com- 
paring the local magnetization below h c i [Fig. [5fc)] and 
in the plateau [Fig. [5jd)] , it is found that the novel de- 
creasing feature of magnetization m T found in case (B) is 
also observed in the present case. The field dependence 
of sublattice magnetization is also studied. As illustrated 
in Fig. [5] for J2/I Ji|=0.5 and 1.0, both m T and ms have 
decreasing regions from h c \ to h C 2- Comparing with the 
sublattice magnetization of case (B), we notice that the 
two sublattice magnetizations that have a decreasing re- 
gion from h c i and h C 2 are those coupled by FM interac- 
tions, and the decreasing behavior is only observed below 



the m=| plateau. 



V. TEMPERATURE DEPENDENCE OF 
SUSCEPTIBILITY AND SPECIFIC HEAT 

From the above results, it can be seen that although 
the three cases entirely exhibit FI ground states, the low- 
lying excitations and magnetic properties are rather dis- 



tinct. Thus, in this section, the temperature dependences 
of zero-field magnetic susceptibility and specific heat are 
explored by the TMRG method.— In the following cal- 
culations, the width of the imaginary time slice is taken 
as e=0.1, and the error caused by the Trotter-Suzuki de- 
composition is less than 10~ 3 . During the TMRG itera- 
tions, 120 and 200 states are retained for the evaluation 
of the susceptibility and specific heat, respectively, and 
the temperature is down to /cBr=0.025| Ji| in general. 
The truncation error is less than 10~ 4 in all calculations. 

The temperature dependence of the susceptibility \ 
and susceptibility temperature product yT for the cases 
are shown in Figs. [3Ja)-(c). For case (A), the susceptibil- 
ity, as shown in the inset of Fig. [TJa), diverges as T— >0 
due to the gapless branch uji.k [Fig. [2fa)]. Upon lower- 
ing temperature, \T decreases to a broad minumum at 
a temperature T m in, and then increases to a peak at low 
temperauture Tp ea it. The minimum of \T is an indicative 
of the Fl-like behavior similar to that in the spin-(i, 1) 
mixed-spin chain. With increasing J2I J\i 2min shifts to 
higher temperatures, corresponding to the enhancement 
of the branches W2,fc and 013^ with the increase of cou- 
pling ratios. Meanwhile, xT min increases for J 2 /Ji<l 
and decreases for J2 / J\ > 1 . The maximum of x^min is 
reached at T m j n =1.25Ji when J2/Ji = l. It is also no- 
ticed that the %T curves for different couplings intersect 
at the same temperature 1.25Ji, as shown by the arrow 
in Fig. [7ja). At low temperature, yT does not diverge, 
like that in the spin-(i, 1) mixed-spin chain^, but has a 
sharp peak, which indicates that x diverges equally or 
slower than y as T— *-03i For the low-temperature peak, 
it is unveiled that T pea k moves to higher temperatures 
with the increase of the height for J2I ' J\<1, while it ap- 
proaches lower temperatures with the height decreasing 
for J2 / J\>1. It can be seen that the finite-temperature 
magnetic properties have transition behaviors with the 
change of the couplings at J2/ Ji=l, which was also noted 
in the ground states in Sec. IV(A). 

For case (B), x also goes to infinity as T— >0 owing 
to the gapless branch u)\ t k [the inset of Fig. EJb)]. As 
shown in Fig. [7][b), xT decreases rapidly to a minimum 
at T m in with decreasing temperature, and then increases 
to a peak at lower temperature, which is quite different 
from that in the case (A). With increasing | J2I/J1, both 
T m ; n and x^min enhance. As indicated by the arrow in 
Fig. [Tfb), xT curves for different couplings also intersect 
at a temperature T^1.15Ji. At low temperature, both 
the peak temperature and height of %T increase with 
increasing l^l/Ji- The peak suggests that x diverges 
equally or slower than ^ as T^0. Different from the 
case (A), the variation of xT in the present case with FM 
coupled pendants does not exhibit a transition behavior. 
It can be seen that J2 has a great impact on the low- 
lying excitations as well as the magnetic properties at 
finite temperature. 

Figure [TJc) illustrates the behaviors of xT for case 
(C). Although x also diverges as T— >Q [the inset of Fig. 
E£c)], xT has rather distinct behaviors from cases (B) 
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FIG. 7: (Color online) Temperature dependence of yT for (a) 
Ji>0, J 2 >0; (b) Ji>0, J 2 <0; and (c) Ji<0, J 2 >0. The insets 
show the susceptibility as a function of temperature. 



and (C) with Ji>0. For J 2 /|Ji|=0.2, \F increases to a 
broad maximum with decreasing temperature, and then 
declines. When 3ij\ Ji|>0.5, a minimum of %T emerges, 
and a small peak appears at a lower temperature. For 
Jil\J\ | >1 j X^ 1 decreases to a minimum with declining 
temperature, showing the AFM feature, and then in- 
creases to a small peak, which is similar to that in the 
case (B). The minimum temperature T m i n also increases 
with enhancing J^/\ J\\. The convergence of yT as T^»0 
indicates that x diverges equally or slower than y as 
T— >0 in this case. Compared with the above cases, no 
intersection of %T is observed for the present case with 
Ji<0. It should be noted that the similar behavior of \T 



FIG. 8: (Color online) Temperature dependence of the specific 
heat C for (a) Ji>0, J 2 >0; (b) Ji>0, J 2 <0; and (c) Ji<0, 
J 2 >0. 



has also been observed in the spin-(|, 1) AFM chain with 
FM next-nearest-neighbor coupling^ which has an anal- 
ogous low- lying excitations. As shown in Fig. [Sfa), the 
gapped magnon branch W2,fc has lower energies than the 
gapless branch 0J\,k f° r large wave momenta k. Thus, the 
low- lying excitations are dominated by W2,fe for small J2. 
With increasing J^/\Ji\, u>2,k enhances and uj\ t k gradu- 
ally dominates the low- lying excitations. The branches 
u)i t k and u>2,k become analogous to that of the case (B) 
[Fig. [3(a)] for large J 2 /| J\\, yielding the behaviors of \T 
for J 2 /\ Ji\>l similar to that of the case (B) [Fig. [7(b)] . 

In Figs. |S(a)-(c), the temperature dependences of the 
specific heat for the three cases are shown explicitly. For 
case (A), the specific heat has a prominent double-peak 
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structure. When J 2 /Ji=0.5, the high-temperature peak 
of specific heat is close to the peak temperature of that in 
the spin-(i,l) mixed-spin chain. With further increas- 
ing J 2 , the low-temperature peak shifts to higher tem- 
peratures when J 2 /Ji<l, while it keeps nearly intact for 
J%j J\>\. Meanwhile, the high-temperature peak con- 
tinuously moves to higher temperatures, which might be 
owing to the enhancement of the gapped branch ioi,k and 

The temperature dependence of specific heat for case 
(B) is shown in Fig. Mb). When | J 2 |M=0.5, the spe- 
cific heat has double peaks, and the high-temperature 
peak is also close to the peak temperature of that in the 
spin-(i,l) mixed-spin chain. Compared with the spe- 
cific heat of the case (A) with J 2 /Ji=0.5 [Fig. [5(a)], it 
can be seen that the high temperature behaviors above 
the high-temperature peak of the two cases agree well 
with each other, but the low-temperature behaviors are 
distinct. With increasing |Ja|/Ji for \J%\/ J\<1, the 
low-temperature peak moves to higher temperature side, 
while the high-temperature peak keeps nearly intact. For 
\Ji\l Ji>1-, the double peaks merge into a single peak, 
which moves to higher temperatures slightly with increas- 
ing IJ2I/J1. Analogous to yT, J2 has also an essential 
effect on the behaviors of specific heat. 

For case (C), the specific heat behaves quite differ- 
ently from the above cases. For J 2 /|Ji |=0.5, the specific 
heat shows a single peak instead of double peaks at low 
temperature. With increasing J 2 below J 2 /|Ji| = l, the 
specific heat below the peak temperature keeps nearly 
unchanged, while the part above the peak tempera- 
ture decreases more slowly, as shown in Fig. [Etc). 
For J 2 /|</i|>l, a high-temperature peak emerges, which 
moves to higher temperatures with increasing J 2 / 1 J\ \ . 
Meanwhile, the behaviors of specific heat below the low- 
temperature peak still retains nearly intact. The low- 
temperature peak seems to be insensitive to J 2 , and is 
dominated by J\. 

For a comparison, we also calculated the thermal quan- 
tities by the LSW theory, which give rise to the similar 
behaviors for the three different cases. The results show 
that yT diverges as T— >0, and the specific heat always 
exhibits double peaks. Although the obtained low-lying 
excitations are helpful to understand the thermodynam- 
ics, the quantitative results obtained from the LSW are 
not so good, which are thus not presented here. 



VI. SUMMARY AND DISCUSSION 

In this paper, the low-lying, magnetic and thermody- 
namic properties of the spin-(i, 1) decorated mixed-spin 
chain with spin-1 pendant spins are systematically stud- 
ied for three cases: (A) J\, J 2 >0; (B) Ji>0, J 2 <0; and 
(C) Ji<0, J 2 >0 by jointly using a few different methods. 

By means of the RSRG analysis, the low-energy ef- 
fective Hamiltonians for each cases in strong and weak 
couplings are obtained. It is found that although the ef- 



fective Hamiltonians are different for three cases, their 
magnon excitations from So to Sq ~ 1 are all FM and 
gapless, which agree with that of the spin-(^, 1) mixed- 
spin chain without pendants. The low-energy dispersions 
of the gapless branch near k=0 are dominated by J\ for 
each case, which is confirmed by the LSW results. 

The low-lying excitations and magnetic properties are 
then investigated by the LSW and DMRG methods, re- 
spectively. The magnon spectra are found to consist of a 
gapless and a gapped branches from Sc to Sc — 1, as well 
as a gapped branch from Sq to Sq + 1, which have differ- 
ent features for three cases. For case (C), two low-energy 
branches have a novel intersection. In a magnetic field, 
case (A) has a m—j plateau, while both cases (B) and 
(C) exhibit two plateaux at ra=\ and |. The low-energy 
gap of case (A) increases almost linearly with increasing 
the coupling ratio, while those of cases (B) and (C) in- 
crease and go to saturation for large |J 2 |, which implies 
that the low-energy gap of the case (A) is mainly scaled 
by J 2 , and those of the cases (B) and (C) are scaled by 
Ji for large |J 2 |. The sublattice magnetization of the 
spins coupled by FM interactions for cases (B) and (C) 
are found to decrease in some regions from h c \ to /i c2 
with the increase of the magnetic field, which may be 
attributed to the competition of the AFM and FM inter- 
actions in a magnetic field. 

The zero-field thermodynamics are also explored by 
means of the TMRG method. It is unveiled that although 
y diverges as T— >0, yT has rather different behaviors for 
each cases. For case (A), yT has a broad minimum and 
a peak at low temperature. The curves of yT for differ- 
ent couplings intersect at a common temperature 1.25 Ji, 
and yT has a transition behavior with the couplings at 
JilJ\=\. For case (B), yT has a narrow minimum and 
a sharp peak at low temperature. The curves of yT for 
different couplings also intersect at a common tempera- 
ture, but yT never show a crossing behavior. For case 
(C), yT has a broad peak for </ 2 /|Ji|<l, and exhibits a 
broad minimum and a peak for J 2 /| Ji|>l, showing two 
distinct features with changing the couplings due to the 
intersection of two low-lying excitations. Compared with 
the spin-(i, 1) mixed-spin chain, there is a common fea- 
ture for the three cases that yT converges as T— »0, which 
implies that y diverges equally or slower than ^ as T^O. 

The specific heat for case (A) has double peaks. 
For case (B), the specific heat has double peaks when 
|^2|/^i<l, which merge into a single peak as | J 2 |/Ji>l. 
For case (C), however, the specific heat has a single 
peak when J 2 / 1 J\ | < 1 , while double peaks emerge when 
Jij\J\ I >1< In a wide range of the coupling for case (C), 
the low-temperature peak appears to be insensitive to J 2 , 
which mainly affects the high-temperature behaviors of 
the specific heat. 

Based on the above results, it can be seen that the 
case (A) of the present system preserves some features of 
the spin-(i 1) mixed-spin chain, while the cases (B) and 
(C) exhibit more novel exotic properties that have not 
been observed in the mixed-spin chains. We expect that 
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the magnetic and thermodynamic properties presented 
in this paper could be tested experimentally in future to 
unveil the effects induced by the pendant spins in the 
mixed-spin chains. 
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